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Abstract. 

"Forbidden" Bragg reflections of iron orthoborate FeaBOe were studied 
theoretically and experimentally in the vicinity the iron K-edge. Their energy 
spectra are explained as resulting from the interference of X-rays scattered at two 
inequivalent crystallographic sites occupied by iron ions. This particular structure 
property gives rise to complex azimuthal dependences of the reflection intensities 
in the pre-edge region as they result from the interplay of site specific dipole- 
quadrupole and quadrupole-quadrupole resonant scattering. Also evidenced is an 
anisotropic character of the absorption spectrum. Self-absorption correction to 
the diffraction data, as well as possible contributions of thermal vibrations and 
magnetic order, are discussed. Particular care is given to extracting clean spectra 
from the data, and it is demonstrated that excellent results can be obtained even 
from measurements that appear corrupted by several effects such as poor crystal 
quality and multiple scattering. 
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1. Introduction 

Resonant X-ray Scattering (RXS) has emerged in the last two decades as a popular 
method for studying local properties of crystals, such as magnetic, charge and orbital 
ordering, thermal vibrations etc. [1, 2, 3, 4]. The anisotropy of the tensor of scattering 
in the vicinity of an absorption edge allows for the observation of Bragg reflections 
that are forbidden in the isotropic case by glide- plane and/or screw-axis symmetries. 
Such reflections were first observed by D.H. and L.K. Templeton [5] and theoretically 
described in [6, 7]. Since then, "forbidden" reflections were studied experimentally in 
many crystals. They are also referred to as pure resonant reflections, because they 
appear only very close to absorption edges owing to resonant scattering processes. 
Because of this energy constraint and of their weakness, they can be measured only 
with synchrotron radiation. 

The sensitivity of RXS to local anisotropy arises from multipole electronic 
transitions and is hence described in terms of scattering tensors of various ranks. 
The strongest contributions to pure resonant reflections are provided by the so 
called dipole-dipole scattering. Beyond this, there are numerous physical phenomena 
contributing to the resonant reflections. Higher-order contributions, like dipole- 
quadrupole [8, 9, 10] or quadruple-quadrupole [11, 12] become important when 
the dipole-dipole contribution vanishes because of symmetry restrictions or physical 
selection rules. In magnetic crystals, the magnetic reflections are mainly provided by 
the antisymmetric part of the dipole-dipole contribution [13, 14, 15]. The co-existence 
of both a magnetic structure and local crystal anisotropy can, in principle, give rise 
to so called "combined" forbidden reflections [16], which have not been observed 
so far. The study of RXS in GdB4 in [17] revealed both time-even and time-odd 
(magnetic) dipole-dipole contributions to forbidden reflections dominating at different 
temperatures. The situation becomes even more complicated by the possibility of 
further effects contributing to forbidden reflections. Contributions associated with 
thermal vibrations and point defects were predicted in [18, 19], and thcrmal-motion- 
induced (TMI) reflections were observed in Ge [20, 21, 22], ZnO [23] and GaN [24]. In 
addition, there can be non-resonant scattering arising from magnetic modulations. 
The interference between non-resonant magnetic scattering and resonant dipole- 
quadrupole and quadrupole-quadrupole channels was found in hematite (q;-Fc2 03) 
[25]. Thus, in many cases, the energy spectra of forbidden reflections can be explained 
as the result of the interference between two or more scattering channels revealing this 
way an interesting physics of electronic interactions. Interference phenomena may be 
even more important when the resonant atoms occupy different crystallographic sites. 
Such situation was observed in magnetite Fe304 [26]. There, the 16(&) position of iron 
provides the hOO, h — 4n + 2 reflections corresponding to the dipole-dipole scattering 
with a maximum close to the absorption edge, while the ions in the 8(a) positions 
allow only for an essential dipole-quadrupole term in the pre-edge. Thus the peaks 
corresponding to the different positions are separated in energy. For garnets, it was 
shown that different atomic positions can provide different forbidden reflections [27]. 

So far, in all studied cases the interference between various scattering processes 
or different atomic positions could be disentangled by taking advantage of their 
manifestation at diflferent energies or thanks to crystallographic selection rules. Here 
we report the case of iron orthoborate (FesBOe), in which scattering from two 
inequivalent crystallographic positions interferes, and this following the same selection 
rules and featuring the same tensor components. Using theoretical analysis and careful 
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data treatment, we demonstrate that RXS is now a technique sufficiently mature to 

study such a complicated case, even in the presence of experimental complications 
such as poor crystal quality, self-absorption and multiple-scattering. 

2. Theoretical 

In kinematic diffraction theory, which is usually applied to RXS, a forbidden reflection 
intensity is simply determined by the resonant structure amplitude F{H), defined by 
the resonant scatterer's position and by the Fourier transform of the atomic scattering 
factor /(r). The latter can be decomposed into: 



where /° is the energy independent non-resonant charge scattering factor and /™ 

accounts for non-resonant magnetic scattering. The last two energy dependent terms 
describe the real and imaginary parts of the RXS and obtain particular importance 
at energies close to an absorption edge. 

The tensor character and the tensorial properties of atomic scattering factors 
near absorption edges were considered in many studies both in a Cartesian approach 
[28, 29] and also using spherical tensors [14, 12, 30, 31]. In the present paper we shall 
use the Cartesian formalism. Witin this framework, the resonant scattering can be 
written in form of a superposition of the interactions in different multipole orders 



where the summation over repeated indices is implied. 

The vectors k, k' and e, e' are the wavevectors and polarizations, respectively, 
of the incident and scattered radiations. Below we shall also use conventional <t and 
TT polarization vectors, normal and parallel, respectively, to the scattering plane. The 
multipolar expansion was extended up to the electric quadrupolar term, omitting 
magnetic multipolar and higher order electric multipolar contributions, which can be 
expected extremely weak [30]. Thus, the resonant scattering tensor consists only of 
the second rank dipole-dipole tensor £), the third rank dipole-quadrupole tensor /, 
and the fourth rank quadrupole-quadrupole tensor Q, all of which possessing tensor 
elements that are sensitive to the incident radiation energy. 

The symmetry properties of these tensors, in the Cartesian formalism, are as 
following: i) the dipole-dipole tensor is symmetric {Djk = D^j) for nonmagnetic 
crystals. For magnetic crystals, the magnetic character manifests itself as an 
antisymmetric part; ii) the time-even dipole-quadrupole tensor can have both a 
symmetric and an antisymmetric part with respect to permutation of the polarization 
indices. It was shown that the symmetric part of the time-even dipolc-riuadrupole 
tensor gives rise to forbidden reflections in Ge and ZnO [8, 20, 21, 23] whereas the 
antisymmetric part causes some peculiarities of the forbidden reflections in Fe203 and 
Cr203 [9, 10, 25] . A time-odd third rank tensor can also occur in the case of a magnetic 
crystal; iii) in non-magnetic crystals, the fourth rank quadrupole-quadrupole tensor 
is symmetric, i.e. invariant against the permutation of polarization or wavevector 
indices [28]. Qjknm Qkjnm Qjkmn — Qkjmn Qnmjk- Thcsc Symmetry 
properties reduce the number of "distinct" cartesian components, for a non- magnetic 
crystal, to 6, 18 and 21 for the dipole-dipole, dipole-quadrupole, and quadrupole- 
quadrupole tensors, respectively [32]. However, spherical formalism shows that the 



f = f + if'^ + f' + if' 
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number of independent components are 6, 15 and 15 respectively. The remaining 
distinct Cartesian components are therefore Unearly related, although their linear 
relationship is not trivial [33]. 

Generally, the intensity of a forbidden reflection is dominated by the dipole-dipole 
contribution. However, when this contribution vanishes, e.g. due to crystal symmetry, 
the higher order terms become the essential ones. They are also important in the pre- 
edge region, where the dipole-dipole contribution is weak: quadrupole transitions are 
reduced in energy due to multi-electron effects such as screening. 

Since the tensorial character of the resonant scattering factor causes the reflection 
intensities to obtain azimuthal dependencies it is essential to identify the non- vanishing 
tensor elements in order to understand the azimuthal properties of forbidden reflections 
and their energy dependencies. Exactly this analytical task is completed in the 
following section for the title compound, FesBOe. 

3. Pure resonant reflections in iron orthoborate 

The symmetry of iron orthoborate structure is described by the space group Pnma 
(no 62) [34]. The Fe cations occupy two octahedral positions: a special one, 4(c) with 
point symmetry m, the mirror plane m is normal to [010], and a general one, 8{d), 
with point symmetry 1 (ie. no symmetry restrictions). The resonant structure factor 
F[H) thus sums the contributions of the 4 and 8, respectively, Fe cations of both 
groups: 

F{H) = F^^'=^ {H) + F*^'*) {H). (3) 

Naturally, the ratio between the partial structure factors, F^^"^ and F^^'^\ varies 
when the scattering factors vary with the reflection indices. Moreover, their relative 
phases vary strongly upon crossing an absorption edge. Below, it will be shown that 
the complex interplay between the scattering at the two Fe positions can produce 
constructive or destructive interference, leading to completely different energy spectra 
for different H = hkl. 

The pure resonant reflections discussed in the present paper are similar to the pure 
nuclear reflections of Mossbauer radiation, but the energy spectra are quite different. 
In the case of nuclear scattering the resonant peaks are very sharp (~ 10~* eV) 
while for electronic transitions the typical width of resonance lines is about 1-10 eV. 
Iron orthoborate has indeed been extensively studied with Mossbauer radiation, and 
the interplay between the two iron positions was evidenced [35, 36]. Because the 
Mossbauer nuclear transition in iron is very sensitive to magnetism, it allows for the 
determination of the magnetic structure of iron orthoborate [37, 38]. In the case of 
RXS, we are dealing with K-edge electronic transitions of the El and E2 types, which 
are only weakly sensitive to magnetism. We will therefore develop our tensor analysis 
treating Fe as non-magnetic. The validity of this approximation will be demonstrated 
by the experiment. 

Next, we determine the contributions to the resonant structure factor, being 
allowed by space group symmetry. We shall first analyse the case of the general 
position 8{d), with general coordinates {x,y,z) in the unit cell , and then proceed 
to the special position 4(c) which second case can be easily deduced from the first 
one. In a general position (i.e. point symmetry 1), the symmetric D, / and Q 
tensors can have up to 6, 18 and 21, respectively, distinct Cartesian components 
[32] , some of which vanishing for particular hkl values. In the case of hOO, h = 2n + l 
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reflections, which are the experimentaUy measured reflections, only one, D^^, of the 
six Cartesian components of the dipole-dipole tensor survives and gives rise to the 
following contribution to the structure factor, where the fractional coordinates (x, y, z) 
refer to the crystal axes (Fig. 1): 

Fdd{h = 2n + 1,Q,Q) ^^D^^cos{2TThx) \ I (4) 





Figure 1. Diagram of the experimental geometry and definition of the coordinate 
system. In this Cartesian system {x,y,z), the vectors fc, fc', <t, tt and tt' are 
written as follows: 

fc =1 fc I ( — sin 9, — cos 6 sini/;, cos 9 cos 

fc' =1 fc I (sin 9, — cos 6 sin i/;, cos 9 cos -0) 

<T = (0, cos i/i, sin ip) 

■K = (cos 0, — sin sin tp, sin 9 cos 

tt' = (cos 9, sin sin ip, — sin S cos jp) 

Similarly, the dipole-quadrupole tensor has a structure factor with only 5 non- 
vanishing distinct Cartesian components at the hOO, h ^ 2n + 1 forbidden reflections: 

^xxz^ Iyyz7 Izxx7 ^zyy 

and Izzz- Describing k and fc' with respect to the crystal {x, y, z) 
coordinate system (Fig. 1), we obtain: 

ixz 

Fd5(/i = 2n + 1, 0, 0) = 4 I A; 1 sm{2Tihx) I -ly, I (5) 



lyz 



with 



ixz = -{Ixxz + Izxx) sin 6* (6) 

iyz = {lyyz ~ Izyy) COS 9 sin l/' (7) 

The term I^zz vanishes because = k'^. For the dipole-quadrupole contribution we 
see that there are both a symmetric and an antisymmetric part. 

Finally, one finds that only 4 of the 21 distinct Cartesian components of the 
quadrupole-quadrupole tensor contribute to the hOO, h ~ 2n + 1 reflections: 



F,,(/i = 2n+ 1,0,0) = 2 I fc 1^ cos(27r/ia;) Qxy (8) 






Qxy 


Qxz 


Qxy 








Qxz 
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Qxy = -Qxyyz COS^ sin 2^' (9) 

Qxz = -Qxxxz Sin^ 6 + COS^ e {Qxzyy Sm^ 4> + Qxzzz COS^ V') (10) 



By substituting the polarizations a = a', tt and w' by their Cartesian expressions 
in equations (4), (5) and (8), we conclude that the x-rays are scattered only in the 
rotated channels, i.e. a — > tt and tt — > a. For the former one, which is the case of the 
experimental setup, we obtain the following contributions: 

F^J{h = 2n+ 1,0,0) = 8r>a;^cos(27r/ia;)cos6isin^ (11) 

F^^ih = 2n + 1, 0, 0) = - 2 I fc I {lyy, - I,yy " I^^, " I,^^) 

X sin(27r/ix) sin20sin V' (12) 
Fq7(/i = 2n + 1, 0, 0) = 2 I fc p cos(27r/ix) 

x[Qi{e)smi; + Q3{e)sm3i;] (13) 

where Qi{9) and QsiO) are linear combinations of Qxxxz, Qxyyz, Qxzzz and Qxzyy 
with coefficients depending on 6 only. This analysis shows that the dipolc-dipolc 
and dipole-quadrupole amplitudes of the scattered radiation possess a simple sin?/' 
azimuthal dependence, while the quadrupole-quadrupole contribution mixes sintp and 
sin 3V' terms. 

The above expressions were derived in the case of the general 8(d) position, for 
which no symmetry applies to the atomic site. We end up with six contributions to 
the forbidden reflections hOQ,h = 2n + 1 : one for the dipolc-dipolc (D^z), one for 
the dipole-quadrupole (a fixed combination of the Ijkn), and four for the quadrupole- 
quadrupole. Let us now analyze the case of ions in 4(c) positions with point symmetry 
m: they lie on the mirror planes {y = ■j and equivalent planes) which transform (xyz) 
into {x, \ —y, z), etc. Thus, all the Cartesian tensor components with indices containing 
odd powers of y vanish. However, equations (11), (12) and (13) do not involve such 
components for the hOO, h — 2n + l reflections so that the scattering description for the 
cations in the 4(c) positions needs as many components as do the cations in general 
8{d) positions. Consequently, for reflections hOO, h = 2n + 1, the resonant amplitude 
scattered by each iron site is described by six parameters, i.e. twelve in total, all of 
which varying with energy. In principle, it would be possible to determine them all, 
by measuring all 6 hOO, h = 2n+l forbidden reflections accessible at the iron K-edge 
energy, at two different azimuths. But it would be extremely tedious, if experimentally 
possible, and it is not the purpose of the present article. Here we want to highlight 
the interference between non-equivalent crystallographic sites. 

The scattering contributions interfere with weights cos(27rft,.T) for the dipolc-dipolc 
and quadrupole-quadrupole, and sin(27r/ix) for the dipole-quadrupole, respectively, 
where x = Xd = 0.41246 for the iron in the 8(d) and x = Xc = 0.12835 for the other 
iron in the 4(c) position [39]. The values of the respective coefficients are given in 
Table 1. Inspecting these values we can see that the dipole-dipole contribution to 
the 300 reflection is mainly provided by Fe in 8(d). We can further see that the 
dipole-dipole contributions from 8{d) and 4(c) to the 300 and 500 reflections have 
opposite signs while those to reflection 700 show comparable values and equal signs. 
Thus, the interference between the radiation scattered by the two iron sites is of the 
same type for the reflections 300 and 500, and opposite for the 700 reflection so that 
signiflcant differences between the energy spectra of the 300, 500 and 700 reflections 
can be expected. 
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{h,k,l) 


cos(27r /iXc) 


sin(27r/ia;c) 


cos(27r/iXd) 


sin(27r/iXd) 


300 


0.0792 


0.9968 


-0.750 


0.661 


500 


0.924 


0.381 


-0.628 


-0.777 


700 


0.7591 


-0.650 


0.803 


-0.595 



Table 1. Weights to the contributions in the hOO, h = 2n + l forbidden reflections 
of the two sites with positions = 0.12835 and Xc = 0.41246. Positions from 
[39]. 



4. Experimental 

The resonant reflections 300, 500 and 700 were measured at the XMaS beamUne 
of the European Synchrotron Radiation Facility (ESRF) at room temperature. 
Complementary measurements of the fluorescence yield and of the resonant reflection 
700 at ambient as well as at low and high temperatures were done at the beamline 116 of 
Diamond Light Source (DLS). The natural linear polarization of the radiation obtained 
at the bending magnet (XMaS) or undulator (116) was used and the energy was tuned 
to the Fe Jf-edge with a Si(lll) monochromator. At both beamlines, the scattered 
radiation was measured with silicon-drift detectors, allowing simultaneous records of 
the elastic scattering and the fluorescence yield in independent channels. Since, as 
shown above, the resonant scattering involves a full u — *■ tt polarization change it was 
not necessary to use a secondary polarization analyzer. For all measurements, the 
same sample was used: a platelet with (±100) faces of a few square millimeters size 
and a thickness of a few hundreds microns. 

In the kinematic diffraction theory, which is usually used for the treatment of 
resonant diffraction, and in case of a sample larger than the beam and much thicker 
than the absorption length of the crystal material, the intensity of the Bragg reflection 
corresponding to the reciprocal vector H = hkl is equal to [40]: 

^ \F{E,e,e'M)\' 
^ ' fi{E,e)+fiiE,e')xg ^ ' 

where the linear absorption coefflcients, /x(£^, e) and n{E, e'), depend on both energy 
E and the polarizations e, e' of the beams, and g is a geometrical factor defined 
as: g = ff^, ct and /3 arc the incident and exit angles with respect to the sample 
surface. In the case of /lOO reflections and a (100) sample surface, g = I. Equation 
(14) assumes that the photon polarization is modified only by the scattering process, 
but not by absorption. Thus, in order to extract \F{E, e, e' , H)\'^ from the RXS 
measurements, the measured intensity must be corrected for self-absorption, which 
requires the determination of the energy dependent absorption and an analysis of its 
anisotropy. 



4-1- Study of the absorption spectrum 

FesBOe has a crystal structure with point group mm,m. The resonant absorption 
is therefore not isotropic. The dipolar component is trichroic (3 independent 
spectra) and the quadrupolar one contains 6 independent spectra. The mixed dipole- 
quadrupole component on the other hand vanishes because the point group permits 
only parity-even absorption events. 

Let us consider the azimuthal dependence of the absorption spectra about 
the [100] direction, with azimuthal reference [001], for any incidence angle 9 and 



Interplay of inequivalent atomic positions in resonant x-ray diffraction 



8 



any polarization. Though the polarization vector is generally not an eigenstate of 
the optical system, wc can make this approximation if the anisotropy is not too 
strong compared to the isotropic absorption. Within this approximation, the linear 
absorption coefficient /U is proportional to the absorption cross-section. According to 
Brouder's formalism [30], fj, can be written as: 

H{E,e,^)=Hnres + miE,9,^) (15) 

{E, 9) + At2(£;, 9) cos(2V') + iH{E. 9) cos(4V') 

where the non-resonant part of the absorption Unres is almost constant across the 
i^-edge, whereas the resonant part at the i^-edgc, ^k, is different for different 
eigcnpolarizations. Since /io and /i2 arc largely ciominatcd by the dipolc contribution 
they can be considered as independent of 9 in case of (t polarized light. The term /X4 
is of quadrupolar origin and occurs only in the pre-cdgc region; ^2 and /X4 are later 
referred to as 2-fold and 4-fold parts, respectively. 

A full energy-azimuth map of fluorescence was measured from the same sample at 
the DLS beamline 116, with the detector at 90° from the incident beam, the incident 
angle at 20° and the polarization pcrpcndicidar to the [100] axis: for simplicity, wc 
shall now refer to this polarization as a, like for the scattering case (Fig. 2, top left). 

For a thick sample, the measured fluorescence intensity is [42]: 

Wl 0, /-^(^) (16) 

where the angular dependence is implicit. Ef = 6.4 keV is the average energy 
of the Kai and Ka2 fluorescence lines. The energy-azimuth map was fitted with 
/io(-B), fJ-iiE), ha{E) and a scale factor as free parameters. Hnres{E) and n{Ef) 
were calculated from the atomic data tables [41]. The final fit parameters and the 
resulting model map are shown in Fig. 2. The measured and calculated maps do not 
show any significative difference, attesting for the good quality of the fit. The 2-fold 
and 4-fold parts of the absorption are clearly evidenced and are respectively 1 and 
2 orders of magnitude lower than the isotropic part. Note that these independent 
spectra are still dependent of 9, and measurements at different incidence angle would 
be needed to evidence it. However, only the quadrupolar part of the absorption 
is 6'-dependent with this polarization, and it is expected to vanish except around the 
pre-edge feature. Ho{E) and ijl2{E) are therefore a good approximation of the isotropic 
and 2-fold components of the dipole absorption for cr-polarized light. 

A similar procedure was applied to the fluorescence data measured along with the 
RXS data at XMaS: each energy-azimuth map was separately fitted. Since, however, 
the data sets were rather noisy, the 2- and 4-fold components could not be extracted. 
The isotropic maps are shown in Fig. 2. Although the forbidden reflections change 
the light polarization from a to tt, as shown in the theoretical section, the diffracted 
intensity is so small compared to the incident intensity that one can infer in first 
approximation that the fluorescence results only from absorption events involving a 
polarized photons. Thus, the case described in the previous paragraph applies, and 
the dipolar absorption should be the same for all /lOO reflections. 

4-2. RXS data reduction 

For the extraction of \F{H)\'^ from the measured data several effects must be taken 
into account. 
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Figure 2. Upper left: Energy-azimuth map of the fluorescence of FeaBOe at 
the Fe K-edge and about the [100] direction. The detector is at 90°, vertical 
to the incident beam, the incident angle is around 20° and the polarization is 
horizontal. Upper right: calculated map from the fit. Lower left: Fit parameters. 
Lower right: Isotropic part of the absorption spectrum, as obtained from fitting 
the fluorescence data. The hkl indices refer to the RXS reflections at which the 
fluorescence data were recorded. Except for reflection 300 the mutual agreement 
of the fitted spectra is reasonably good. The 300 data set was rather poor due 
to the combination of the low incidence angle, the large beam size and a sample 
inhomogeneity. Except for the pre-edge region the general shape of the absorption 
coefficient is in agreement with the FDMNES calculations (see Section 5). For 
comparison, also the values calculated from Henke's tables [41] are shown. 

The first one is absorption effect wfiicli distorts tfie energy spectra. In principie, 
tlie absorption fias also an angular dependence, but as shown in the previous section 
ignoring it is not a bad approximation. Hence, an effective absorption spectrum was 
calculated by averaging the 4 reasonably isotropic absorption curves derived from the 
fluorescence measured at the reflections 500, 700, 710 and 710 (300 was not considered 
for the reasons given above). According to Equation (14), the measured data were 
multiplied by IJ.{E, e) + guiE, e') w (1 + g){finres + IJ-o{E)), where g is 1 for the hOO 
reflections. Of course, correcting the data with an isotropic absorption spectrum and 
neglecting its azimuthal dependence introduces errors. From Figure 2, bottom left 
panel, we estimate that these can locally amount to 15% (particularly in the pre-edge 
region), and to 5% at energies above the edge. 

Another problem is the inhomogeneity of the sample, which leads to strong 
variations of the rocking curve profiles obtained in azimuthal scans: the ratio 
of the maximum height of the rocking curve to its integrated intensity is not 
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Figure 3. Left: rocking curves of tlie reflection 500 at two different azimutfis, 
illustrating the variability of the shape (hence the crystal quality). Right: 
maximum intensity to integrated intensity ratio of the 500 reflection rocking curves 
versus azimuth. The intensities recorded at 4 different energies show that this ratio 
is essentially energy independent. 



constant. Since the energy-azimuth maps were recorded "on the fly", i. e. without 
integrating each rocking curve, we had to correct for the azimuthal dependence of the 
maximum/integrated intensity ratio of rocking curves, that were recorded at a few 
sample energies (Fig. 3). As expected, these ratios seem to be energy independent 
within the small energy range of the spectra. 

The last effect to correct for is multiple scattering. In fact, most of the energy 
spectra are severely contaminated by multiple scattering (Fig 4, top left). To correct 
for this, we eliminated all experimental points whose intensities relative to those of 
their azimuthal neighbors exceeded a threshold value which was chosen as the standard 
deviation of the next-neighbour intensity variation. The eliminated points were then 
replaced by interpolations from the remaining points. This method takes advantage of 
the small azimuthal step size (5°) of the experimental data and most of the multiple 
scattering could be removed, the remaining contributions being too small to have 
significant impact on subsequent data analysis (Fig 4, top right). 

The resulting two-dimensional azimuth-energy maps of intensities of the forbidden 
reflections 300, 500 and 700 are shown in Figure 4. 

In a further step, according to the theoretical predictions (Eq. (11), (12) and 
(13)), the corrected data I{E,ip) were fitted with the model function: 

liE, V) = a{E) sin + b{E)e"^'^^^ sin(3V') ^ (17) 

where a{E), b{E) and $(£') are real positive parameters (Fig. 5). The differences 
between the 'corrected' maps and the fitted maps are essentially due to the residual 
multiple scattering. As expected, the fit curves show a strong 2-fold component, 
mainly due to the dipole-dipole scattering which largely hides the 2-fold dipole- 
quadrupole scattering as well as the 2-fold part of the quadrupole-quadrupole 
scattering. The 4-fold component of the quadrupole-quadrupole scattering appears 
only in the pre-edge region as expected. The small features at higher energies may be 
explained by the isotropic absorption correction: the 2-fold aspects of the scattering 
(|sin?/;|^) and of the absorption (cos 2^) yield a weak 4-fold contribution, like the 
interference between the sin^ and sin3'0 terms in the scattering amplitude. This 
contribution appears on the main peaks of the reflections 300 and 500, where both 



Interplay of inequivalent atomic positions in resonant x-ray diffraction 11 



300 300 




-150 -100 -50 50 100 150 -150 -100 -50 50 100 150 

Azimulh (deg) Azimuth (deg) 



Figure 4. Energy-azimuth maps of the 300 (top), 500 (bottom left) and 700 
(bottom right) forbidden reflections. For the 300, both raw data (left) and 
corrected data (right) are shown. Note the distortion of the azimuthal dependence 
by the poor crystal quality, the distortion of the energy dependence by the self- 
absorption, and the dense multiple scattering. For the 500 and 700, only corrected 
data are shown. 

intensity and absorption are large. Elsewhere it is within the noise of the measurement. 
Its weakness proves that our approximation of isotropic absorption sufSces for the 
correction of the RXS data. The phase parameter ^{E) remains close to 0[7r], except 
for reflection 700 for which the 4- fold feature seems to be shifted in phase by 7r/2[7r] 
against the 2-fold one. However, the data for the fit are too noisy to state this with 
confidence. 

4-3. Thermal- Motion- Induced scattering and magnetic scattering 

So far, our analysis is based on the hypothesis that the forbidden reflections are 
exclusively due to time-even scattering of the atoms in their nominal positions, i.e. 
we neglect contributions of Thermal Motion Induced Scattering (TMI). TMI is the 
derivative of the scattering tensor of the atom at rest with respect to its position. Due 
to the predominating dipole-dipole scattering, its main contribution is therefore a third 
rank tensor with the same symmetry constraints as those of the symmetric dipole- 
quadrupole scattering tensor, and, being related to thermal displacements of atoms. 
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Figure 5. Top: selected azimuthal dependences and their fits against equation 
(17). A liint of asymmetry, visible in the pre-edge of reflections 300 and 700, 
can be attributed to non-resonant magnetic scattering. Bottom: fit parameters 
a{E) and b{E) of the energy-azimuth maps of Figure 4 with the model given in 
equation (17). Left: 300, middle: 500, right: 700. 



TMI almost vanishes at low temperatures. Therefore, we measured the reflection 700 
also at low temperature. The energy spectra at room temperature (300 K) and at low 
temperature (8 K) are almost identical (Fig. 6, left). This finding proves that TMI 
must be small compared to the tensor of scattering at rest. 




Figure 6. Energy spectra of the forbidden reflection 700 at low temperature 
(8 K), room temperature (300 K) and elevated temperature (550 K). The 
normalized intensities without absorption corrections are presented. Note that 
the data at room temperature were taken at XMaS, whereas the data 8 and 
550 K were measured at DLS 116. Left: Azimuth ip = 90°. Right: Azimuth 
^ = -72° = -90° + 18° at 550 K and i/i = -108° = -90° - 18° at 300 K. 

Also SO far, we have considered the crystal as non magnetic although it is 
weakly ferromagnetic with a paramagnetic transition at Tc = 508 K. Despite of 
the small ferromagneticity, most of its ordered moments follow an antiferromagnetic 
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order with propagation vector [000]: the spins form two ferromagnetic sub-layers in 
the (a.6)-plane, and the coupling between the two sub-layers, along the c-axis, is 
antiferromagnetic [43]. It undergoes a spin reorientation transition at Tsr = 415 
K, when the antiferromagnetic and ferromagnetic axes swap their directions: the 
antiferromagnetic axis goes from [001] below Tsb. to [100] above Tsb. [44]. As a 
consequence of the antiferromagnetic component, the magnetic structure factor does 
not vanish for the forbidden reflections hOO, h = 2n + 1, i.e. both resonant and non- 
resonant magnetic scattering can contribute to their intensities. 

The magnetic moments are carried by the Fe'^+ cations, which, in iron oxides, 
are normally in high spin state (5 = 5/2 and L = 0) [45]. Since resonant electric 
transitions are only sensitive to orbital moments (and to the spin moments via 
spin-orbit coupling), a resonant magnetic scattering contribution to the forbidden 
reflections can be expected very weak, but a contribution of non-resonant magnetic 
scattering cannot be ruled out. Indeed, the experimental spectra derived at room 
temperature and low temperature do not vanish below the absorption edge, which 
suggests measurable effects of non-resonant magnetic scattering (Figs. 5 and 6). This 
assumption of non-resonant magnetic scattering is confirmed by the measurements 
performed at 550 K, i.e. well above the paramagnetic transition, which failed to yield 
any measurable signals below the edge (Fig. 6, right). Close to and above the edge, 
the strong absorption reduces the scattering volume in such way that non-resonant 
contributions become small compared to the resonance enhanced contributions. Since 
our aim is to observe and demonstrate the interferences of the RXS from the two 
iron sites, we will not further focus on the non-resonant region and ignore magnetic 
scattering in the ensuing data analysis. 

5. Fitting the energy and azimuthal dependencies 

Our symmetry analysis has shown that the structure factors of the hOO, h = 2n + 1 
reflections can be conveniently described by 12 tensor components (6 for each site), 
but azimuthal scans allow only for extracting 2 independent spectra per reflection. 
Therefore, the experimental results can only be interpreted by numerical simulations. 
To this purpose, the FDMNES code [46] was used for calculating the absorption 
cross-sections and the structure factors \F{H)\'^ of the forbidden reflections. This 
information served to yield data that can be compared with the experimental 
intensities, after correction. Using the multiple scattering option of the code, with 
the Hedin-Lundquist potential, it was attempted to optimize the model parameters 
with respect to best agreements with the experimental values of I{H) and ij,{E). 

Let us refer to the strong reflections above 7.120 keV as main peaks and to the 
weak peaks below 7.120 keV as pre-edge peaks. The numerical simulation of the 
various contributions to the tensor atomic factors show that while the main peaks 
are mainly due to the dipole-dipole contributions from both of the iron cations, the 
pre-edge peaks can be described by higher order terms producing complex azimuthal 
and energy structures. Despite of the dipole-dipole origin of the main peaks, their 
intensity modelling is not simple, just because of the interference of the radiation 
scattered at the two Fe positions. It follows from Table 1, that 4(c) and 8(d) structure 
factors interfere in very different ways for the 300, 500 and 700 reflections. Indeed, 
the experimentally observed 300 and 500 reflections are strong, contrary to the weak 
700, and the 700 energy spectrum looks rather different from those of 300 and 500. 

In the course of simulations, the main peaks of the 3 measured reflection intensities 
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Figure 7. Energy spectra of the reflections 300, 500 and 700 at the azimuth 
j/j = 90°: FDMNES calculations versus experimental data. 



were fitted first. Neglecting higfier order terms in their energy range, the model 
intensities depend on two parameters only: D^f ^ and dH^K Table 1 shows that the 
dipole-dipole component of the 300 intensity is mainly provided by the 8{d) sites, 
while the contribution of the 4(c) sites is very small. This circumstance allows for 
a convenient extraction of Dx^f \ whereupon ot^^'' could be extracted from the 500 
spectrum, once Dx^f^ was known. At last, the 700, in which the radiation interference 
between the iron sites is mainly destructive, was used for fine tuning, because the 
destructive interference makes modelling extremely sensitive to the ratio of both 
parameters. A satisfactory description of the 700 spectrum could only be achieved 
upon introducing an energy shift of 0.7 eV between both contributions. This chemical 
shift can be rationalized by the different environments of the iron cations occupying 
two inequivalent sites. Then, it was checked that the same set of parameters allows 
for reasonably good 'fits' of all 3 reflections (Fig. 7) as well as of the absorption (Fig. 
2). 

The key step of the simulation procedure is the convolution that describes the 
spectral broadening r(£^), due to the core- and valence electron lifetimes in the excited 
state. This function being unknown, FDMNES supposes T{E) to grow smoothly, e.g. 
Hke 

r = Tiiole + I 



11 /tt r„, / 1 

— I — arctan — 



(18) 



2 TT y 3 Elarge 

where e = ^-^F.r-rai _ ^j^g 5 parameters Thole, T^, Ecent, Epermi: and Eiarge, 
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that provided the least differences from the experimental data, were obtained with 
the optimization procedure of the FDMNES code, yielding: Ep^rmi — — 6.5ey, 
r„, = E,,^t = 14ey, Eiarge = lOeV, Vhoie = 2el^. 

Since the main peaks, near and above the edge, result from the dipole-dipole 
scattering their intensities observed upon rotation about the scattering vector display 
a stable 2-fold azimuthal dependence. In the pre-edge region, however, the features 
observed result from the interplay between the dipole-quadrupole and quadrupole- 
quadrupole terms so that more complex azimuthal dependencies are displayed which, 
in addition, vary strongly with energy (Fig. 8). Their intensity values are too small 
to allow for detailed modelling, but the FDMNES calculations using the previously 
determined convolution parameters reproduce the general shape of the data reasonably 
well (Fig. 8). 



300, E=7.113 keV x 10"^ 300, E=7.1 1 35 keV 




Figure 8. Azimuthal dependence of the reflection 300 at 7.113 keV, 7.1135 
keV, 7.114 keV and 7.12 keV. Comparison between experimental data ('data'), 
fit against equation (17) ('fit'), and simulation with FDMNES ('FDMNES'). The 
simulation fulfills equation 17, but the coefficients are not the best fit to the data. 

Figure 9 shows the calculated dipole-quadrupole and the quadrupole-quadrupole 
contributions to reflection 300 of each iron position, at 7.114 keV and 7.115 keV. 
Although these energies seem very close to each other, the various contributions 
show strong energy dependent changes of azimuthal behaviour, which affect the total 
intensity. 
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Figure 9. Azimuthal dependencies of tlie dipole-quadrupole and quadrupole- 
quadrupole contributions to the 300 reflection at 7.114 keV (left) and 7.115 
keV (riglit). The FDMNES convolution parameters are the same as previously 
determined. 



6. Conclusions 

Resonant X-ray diffraction near the iron K-edge was studied in iron ortlioborate, botli 
theoretically and experimentally. The complex energy and azimuthal dependencies of 
the 300, 500 and 700 forbidden reflection intensities are explained as resulting from 
the interference between the radiation scattered from the two crystallographically 
inequivalent iron positions: 8(d) and 4(c). In the dipole-dipole approximation, the 
300 intensity results mainly from Fe at the 8{d) position, whereas the 500 and 
700 intensities result respectively from constructive and destructive interferences of 
radiation scattered from Fe at both crystallographic positions. The features in the pre- 
edge region can be explained by the interplay of dipole-quadrupole and quadrupole- 
quadrupole processes. Due to the low point symmetry of the iron sites, the large 
number of parameters required for the full description of the resonant scattering 
prevented their complete determination. Nevertheless, the dipole-dipole components 
could be extracted from the large part of the energy spectra where higher rank 
processes are comparatively weak, i.e. near and above the if -edge. Theoretical 
calculations with the FDMNES code and subsequent fits of the experimental spectra 
suggest an energy difference of 0.7 eV between the K-edges of the two non-equivalent 
iron atoms. This chemical shift is attributed to the different crystallographic 
environments, to which RXS is very sensitive. The measurements of the 700 reflection, 
whose net intensity and energy spectrum were found stable over a wide temperature 
range, indicate the absence of significant contributions of magnetic and/or thermal- 
motion-induced scattering. 

The investigations demonstrate also the remarkable sensitivity of RXS to the 
interplay of various radiation contributions, and simultaneously that a thorough 
analysis is possible despite its apparent non-triviality. However, only studies of 
both energy and azimuthal dependencies of several forbidden reflection intensities 
provide the possibility to distinguish contributions from different crystallographic 
sites and hence to study differences in the electronic interactions of atoms thereon. 
Finally, it should be noted that excellent experimental data could be extracted 
from measurements that were seriously affected by the poor crystal quality, frequent 
multiple scattering and strong self absorption. This finding demonstrates as in earlier 
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studies, e.g. [47], that RXS as a local probe is neither restricted to using nearly perfect 
crystals nor to especially favourable diffraction geometries. 
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